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In models with extended scalar sectors with several Higgs multiplets, such as Higgs triplet models, 
the Z gauge boson can have nondiagonal couplings to charged Higgs bosons. In a model-independent 
way, we study the potential contributions arising from such theories to the CP-vioIating trilinear 
neutral gauge boson couplings ZZ*"f, ZZ'f, and ZZZ*, which are parametrized by four form 
factors: /ip2, /4 , and ff, respectively. Such form factors can only be induced if there are at least 
two nondegenerate charged Higgs scalars and an imaginary phase in the coupling constants. For 
the masses of the charged scalar bosons we consider values above 300 GeV and find that the form 
factors can reach the following orders of magnitude: |/if| ~ 10~® — 10“'*, |/if| ~ 10“^ — 10“®, 

1/4 I ~ 10“® — 10“^ and |//| ~ 10~® — 10“®, though there could be an additional suppression factor 
arising from the coupling constants. We also find that the form factors decouple at high energies 
and are not very sensitive to a change in the masses of the charged scalar bosons. Apart from a 
proportionality factor, our results for the form factor, associated with the ZZZ* vertex, is of 
the same order of magnitude than that induced via nondiagonal neutral scalar boson couplings in 
the framework of a two-Higgs doublet model. 

PACS numbers: 12.60.Fr,14.70.Hp,14.80.Fd 


I. INTRODUCTION 

Pair production of electroweak gauge bosons {W, Z and 7) provides an import test for the gauge sector of the 
standard model (SM), and opens up the possibility for observing new physics phenomena associated with such particles. 
In this context the production of a pair of neutral gauge bosons, or ZZ, was studied at LEP 0,0 , Tevatron M 
and the LHC p-0. The main production mechanism proceeds at leading order via the t and u channels, but regardless 
of the production process, the experimental measurements were found to be consistent with the SM predictions 01 - 
Any deviation in the corresponding cross sections could be a hint of new physics effects, such as new heavy particles 
[lol| . new interactions [HI,[13, etc. In particular, Zj and ZZ production could allow us to study the trilinear neutral 
gauge boson couplings (TNGBCs) ZViVj, whose framework is best discussed via the effective Lagrangian approach, 
as in Ref. [l^l, where the lowest-dimension effective operators inducing the off-shell TNGBcs were presented. Since 
the 777 coupling is forbidden by Furry’s theorem, experimentalists have focused their attention on the .^77, ZZ^ and 
ZZZ couplings. On the other hand, Landau-Yang’s theorem forbids any TNGBC with three on-shell gauge bosons, 
so there are only two vertex functions describing four distinct TNGBCs with one off-shell gauge boson: ZV*j and 
ZZV* {V = 7, Z) [13 • The most general ZV*j vertex function fulfilling both Lorentz and electromagnetic gauge 
invariance is parametrized in terms of four forms factors hY (* = I, 2, 3,4), whereas the ZZV* vertex function is 
parametrized by only two forms factors /T (j = 4,5). While fY are CP violating, and fY are CP 

conserving. Any TNGBC is zero at the tree level in the SM or any renormalizable extension, only the CP-conserving 
ones are non-vanishing at the one-loop level in the SM via the fermion triangle [l5l| - whose contribution is highly 
suppressed even in the presence of a fourth fermion family [l^ . The respective experimental bounds on TNGBCs are 
expected to achieve a significant improvement in the LHC in the forthcoming years. Thus, it is worth studying any 
possible contribution to these couplings at a high energy scale A via effective Lagrangian. 

Although the CP violation observed in the K meson system can be explained by the Cabibbo-Kobayashi-Maskawa 
matrix complex phase, the SM does not predicts enough CP violating effects to explain the current matter-antimatter 
asymmetry in the universe. Consequently, other sources of CP violation are necessary. Such CP-violating effects 
may show up via TNGBCs, which are thus worth studying. Along these lines, it was shown that the radiative decay 
Z —/r+/r“7 may be sensitive to both CP-conserving and CP-violating ZV*^ and ZZV* TNGBCs. Furthermore, such 
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a process may also be useful to put constraints on the CP-violating forms factors in future linear collider experiments 
(l7l| . The current experimental bounds on TNGBCs reported by the PDG collaboration [T^ come from a combination 
of LEP measurements 0 , where the experimental results and the individual analyses are based in reports between 
1999 and 2001. However, the L3 collaboration has updated their analyses, resulting in the most restrictive bound on 
up to date [l|. On the other hand, a few results have been reported by the GMS ii and ATLAS collaborations 
The current most stringent bound on fY was obtained by CMS Q, based on data collected in 2010 and 2011 
at y/s = 7 TeV with an integrated luminosity of 5.0 ±0.1 fb“^. A more accurate analysis is expected by both the 
ATLAS and the CMS collaborations in the forthcoming years. The most stringent bounds on the TNGBCs form 
factors are shown in Table HI 


Experiment 

Limit 

L3 

L3 

CMS [8j 
CMS [8] 

-0.153 < /if < 0.141 
-0.087 < /if < 0.079 
-0.011 < fi < 0.012 
-0.013 < f2 < 0.015 


TABLE I: The current most stringent limits on CP-violating TNGBCs. 


Contributions to TNGBCs have been previously studied in the context of the minimal supersymmetric standard 
model (MSSM) and the littlest Higgs model j2l|, focusing only on the CP-conserving form factors. Moreover, 
the CP-violating ZZZ* form factors were studied in the framework of two-Higgs doublet model (THDM), where the 
respective contributions are induced via nondiagonal complex couplings arising in the neutral scalar sector 22|. In 


the framework of several extensions of the SM, flavor change is allowed in the fermion sector via tree level neutral 
currents. If the respective coupling constants contain an imaginary phase, they can induce CP violating TNGBCs at 
the one-loop level. Another possibility arises when tree level nondiagonal complex couplings Z^f^J appear in the 
scalar sector, with charged scalar bosons. In such a case, non-degenerate charged scalar bosons and a complex 
mixing matrix are necessary to induce TNGBCs. In this work we will present an analysis on the CP-violating TNGBCs 
ZZ*j, ZZj* and ZZZ*. In particular we will focus on the one-loop level contributions from nondiagonal complex 
couplings arising in the charged scalar sector of a SM extension. Our analysis will be rather general as we will use the 
effective Lagrangian approach. We have organized our presentation as follows. The vertex functions for the TNGBCs 
and the effective Lagrangian from which they arise are shown in Sec. H. Section HI is devoted to the analytical results 
for the one-loop calculation. Numerical results and discussion are presented in Sec. IV. Finally, the conclusions and 
outlook are presented in Sec. V. 


II. TRILINEAR NEUTRAL GAUGE BOSON COUPLINGS 


The most general effective Lagrangian describing TNGBCs ZViVj contains both CP-even and CP-odd terms 


Furthermore, if all of the gauge bosons are taken off-shell, there are both scalar and transverse structures 
However, the scalar terms vanish when on-shell conditions are considered for the ZV*^ and ZZV* couplings 
Thus, the effective Lagrangians describing such couplings can be written as [2J: 


Czv^^ = -W - [hl{d<^F^Y + -\(hl{do.dpdPFpY ± ± ml)ZY\)z^F>^^ 

TTt ^ V TTh ^ \ / 

-[hlidpFP^) + hl{dpZP>^)]Z‘^F^^ + + hf[{d‘^ + ml)dPz^Y)ZpF^c], ( 1 ) 

Fzzv^ = ^{- + fi{d^Z>*^)]Z^id‘^Zp) + [fUd'^F^Y + . (2) 

Here 1^^ = /2, with = 9^14 — d,yV^ standing for the stress tensor of the neutral gauge boson. The 

operators associated with 4 have dimension eight, but the remaining ones are of dimension six. While hY 2 and fY 
are CP-odd, and fY are CP-even. From the effective Lagrangians (HD and 0, we obtain the vertex functions 
fe(pi,P2,(?) for the ZV*j and ZZV* couplings respecting Lorentz covariance, C/em(l) gauge invariance, and 
Bose symmetry, which are given by 00: 
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pa/3/i 


{Pi,P2,q) 


pa/3/i 


.(pi,P2,g) 


, {pI - mlr) 


hXiq^g'^^ - q^g^n + ^P2 {{.q-P2)9^^ - q^pX) 


-hle^<^^Pp2p - ^p^e^>^Py2pqo 


m% 


{flig^g^^ + <i^9n - fXe>^'^^'’ipi -P2)p) , 


( 3 ) 

( 4 ) 


where my is the mass of the off-shell V gauge boson, and the overall factor {p^ — to|) in Eq. ([3]) is a consequence of 
Bose symmetry. The four-momenta of the gauge bosons are defined in Figure [1] 


V/{q) 



FIG. 1: Vertex function ie F^^^y. (pi,P 2 , ?) for the ZViVj {Vij = Z, 7 ) couplings. All the four-momenta are outgoing. 

In the next section we will consider a SM extension including complex nondiagonal couplings in the charged scalar 
sector. We will show that only the CP-odd form factors are induced by this mechanism. 


III. ZZ*7 AND ZZV* CP-ODD FORM FACTORS 

We consider a renormalizable effective theory with several physical charged scalars bosons.^ We will assume that 
there are trilinear vertices which can induce TNGBCs at the one-loop level via triangle diagrams. The 

Lagrangian describing such trilinear interactions between the neutral gauge bosons and the charged scalar bosons can 
be written as C = Cd + Cnd, where Cd describes the diagonal couplings and Cnd the nondiagonal ones. The gauge 
invariant structure of these Lagrangians is as follows 

i i 

Cnd = -FH.c. 

where the coupling constants gf^ are real due to the Lagrangian hermicity, whereas could be complex. From 
here we extract the Feynman rules for the interaction between neutral gauge bosons and the charged scalars bosons, 
{ki)^~{k 2 ), which has the generic form iegYj{ki — where all the four-momenta are outgoing. For F = 7 , 
there are only diagonal couplings and g]^ = is the charge of the scalar boson in units of e. 

It is worth mentioning that in the class of theories we are considering, quartic vertices of the kind 
would also appear. In principle, this class of vertices can contribute to TNGBCs at the one-loop level through bubble 
diagrams such as the one depicted in Figure^ However since the renormalizable quartic vertex would be proportional 
to the metric tensor 5 “^, the amplitude arising from this of diagram can be written as 


( 5 ) 

( 6 ) 


^ For instance the low energy effective theory arising from the 331 model or a Higgs triplet model. 
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aPfj. 

bubble 


d^k i2k + qrg‘-0 ^ 

(27r)'° (fc2 - mf)((fc + g)2 - to|) ^ ^ 


(7) 


Therefore this diagram does not contribute to our TNGBC form factors. All other possible diagrams obtained by 
permuting the Z gauge bosons contribute with terms that can be dropped when the transversality condition for the 
on-shell gauge bosons are taken into account and also when considering that the virtual gauge boson is attached to a 
conserved current. 



Z“(pi) Zf^{p2) 

WVWVWVAAAAAAAAAA 

FIG. 2: Illustrative bubble diagram arising from the quartic vertex which can appear in the class of theories we are 

interested in. Since this vertex is proportional to the metric tensor g°‘^, the amplitude of this diagram is proportional to q'' 
and does not contribute to our CP-violating TNGBC form factors, which neither receive contributions from all other bubble 
diagrams obtained by exchanging the Z gauge bosons. 


We now present our results for the one-loop contributions to the ZV*^ and ZZV* couplings, where V* stands 
for an off-shell gauge boson. Since only the Z gauge boson has nondiagonal couplings, the .^ 7*7 form 

factor will not be induced at the one-loop level via this mechanism. Also, no CP-even form factor is induced via 
scalar couplings since the Levi-Civitta tensor cannot be generated this way. We will set Q<i> = — 1, but the results 
will also be valid for the contribution of a pair of doubly charged scalar bosons, though in this case there will be 
an additional proportionality factor of 2 appearing in the ZZ*^ and ZZ^* vertices. Once the amplitude for each 
Feynman diagram was written down, the Passarino-Veltman method was applied to obtain scalar integrals, which are 
suitable for numerical evaluation. 


A. ZZ*'y coupling 

The Feynman diagrams contributing to the Za{pi)Zp(p 2 )^^{ci) coupling are shown in Figure[3] It is interesting to 
note that diagram (b) is required by Bose symmetry, whereas diagrams (c) and (d), which involve the exchange of 
the virtual scalars bosons, are necessary to cancel out ultraviolet divergences. After the mass-shell and transversality 
conditions for the gauge bosons are considered, we obtain the following results 






TO|lm . 

67 r 2 (TO 2 -p 2)3 (pI - "*!) “ pI) - m] - m\) Cijjipj) 

- m]) {{ml - pI){B,j(0) - 1 ) + 3{ml + pl)B,j{pl) - 6 TO|By(m|)) - (f O j) (8) 

6,^2(^2 (P^ ~ - im]{ml - pI) {ml -bpl - 2 {mj - to|)) 

]^{m} - ml) (2(m| -p^)(By(0) - 1) -b 3(to| -b ‘ipl)B,j{pl) - 3{3ml + pl)B,j{ml)) - (i ^ j) | 

(9) 


where p 2 is the four-momentum of the off-shell Z gauge boson, mij stand for the masses of the charged scalar 
bosons and the coefficient 2 ilm(g^ 5 ?*) = gfjgf* — gf/gfi contains the imaginary phase that induces CP violation, 
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FIG. 3: One-loop level Feynman diagrams inducing the ZZ*'^ (^ = 7 ) coupling. 


which is necessary to obtain nonzero results and is consistent with the Lorentz structure of the Lagrangian. We 
also have introduced the shorthand notation Bab{<?) = Bo{c^Cabdpl) = m^) and 

Cabciqd = C'o(to|, 9^, Wai ^c)) where Bq and Co stand for Passarino-Veltman scalar functions. From the 
above expressions, it is evident that the form factors vanish when the masses of the charged scalar bosons are 
degenerate. It is also straightforward to show that ultraviolet divergences cancel out. 


B. ZZ^* coupling 

The Feynman diagrams for the Za(pi)Z 0 [p 2 )l*^{(l) couplings are similar to those inducing the ZZ*^ coupling, but 
in this case the photon is off-shell. The corresponding form factor is given by: 




UlF' Im i 

12m|7r2g2(4TO| _ ^2)2 “ "*?) (l2m|Bij (TO|) + (4 to| - “ ^y (O))) 

-I- (rUj -I- 3m| -|- 7 to| — -|- 6 — 2m%{8rrif — 3to|)^ Bu{q^) 

- 6 {mj-rrd^- m|) - 2 m} {m} + m%) + {m} - m?zf + rn}q^^ C\ji{qd - (z O j)|, (10) 

where q is the photon four-momentum. As expected, this form factor is proportional to Im( 5 ^g?*) and vanishes when 

mi = mj. 


C. ZZZ* coupling 


Because of Bose symmetry, in the case of the ZZZ* coupling, there are several more diagrams than those inducing 
the ZZ*^ vertex. We present in Figure U] the generic Feynman diagrams from which all diagrams inducing the 
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^“{Pl) (3) Z'‘{q} 


FIG. 4: Generic one-loop Feynman diagrams for the ZZZ* coupling. The complete set of diagrams is obtained by permuting 
Za{pi) and Zp{p 2 ) and exchanging the charged scalars bosons. 


ZZZ* vertex can be generated. Notice that diagrams (2) and (3) are obtained from diagram (1) after performing 
the permutations Za{pi) O Z^(q), and Zp{p 2 ) O Z^{q), respectively. Additional diagrams are obtained from these 
diagrams following a similar procedure as that described in Fig. [31 for each one of the Feynman diagrams of Fig. 
[4] there are three more diagrams that are obtained similarly as diagrams (b)-(d) of Fig. [31 which are obtained from 
diagram (a) by permuting the Z gauge bosons and exchanging the charged scalars. Therefore, there are a total of 
twelve Feynman diagrams for the ZZZ* coupling. By using the appropriate simplifications, the amplitude of each 
diagram of Figure [4] reduces to those of the ZZ*^ and ZZ^* couplings. The diagrams (2) and (3) of Fig. [4] not only 
are required by Bose symmetry but, once their amplitudes are added up, ultraviolet divergences cancel out. After the 
Passarino-Veltman method is applied, we obtain the following result 


X 

+ 


127r^SwqHq^-ml){q^-4m%)^ 

(rrii + 3m| -b 7m|) -b 6 {ml - mff - - 6 to| - q"^^ Bu{q^) 

{m\q^ (lOm| - I'imi — 3to|) -b 6 m%{mf - m|) {ml - -b 2m|) -b q^ {‘im'f - m|)) Bii{m%) 
(m^ - m"^) {2q^ (4m| - q^) (1 - 5^(0)) - 3q^B,j{q^) - 3m| (4m| - Bij{m\)) 

6 {ml - -b 2 m|) (^m%q'^ (m| - 3mf - m^) + m% ("i? “ ^njiq^) 

Qq^ {mf - TOj - m|) (mf - 2mj (to| -b to|) -b (m| - m|)^ -b C^ji{q^) - (f O j)|, ( 11 ) 


where q is now the four-momentum of the off-shell Z boson. We note that all the properties discussed above are also 
present in this form factor. In the next section we will evaluate the CP-violating TNGBCs for illustrative values of 
the charged scalar boson masses and the four-momentum of the virtual gauge boson. 
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IV. NUMERICAL RESULTS AND DISCUSSION 


While the diagonal couplings can appear in several extensions of the SM at the tree-level, the presence of the 

nondiagonal couplings is less common, but they can be induced indeed within a more general renormalizable 

theory. In order to analyze the CP violating TNGBC form factors, we will not consider specific values for nor 
Therefore, the masses of the charged scalar bosons rriij and the four-momentum of the virtual gauge 
boson will be the only free parameters involved in our analysis. 

A particle with the properties of the SM Higgs boson with a mass of 125 GeV has been finally discovered at 
the LHC 
ATLAS [2313 


26 l| . and the search for singly ($'*") and doubly charged scalar bosons is still underway by the 

and CMS [3l|, 113 collaborations. However, no evidence for the existence of such scalars bosons has 


been found up to date. Based on data collected in 2011 , the ATLAS collaboration performed a model independent 
analysis to search for a light charged Higgs boson with a mass in the range 90-160 GeV [13] ■ Independently, the GMS 
collaboration reported a search for the charged Higgs boson of the MSSM with a mass ranging from 80 to 160 GeV 
[Hlj . This collaboration also reported a lower bound on the doubly charged Higgs boson mass between 204 and 449 
GeV from the processes pp —?> $++$ —>• ^ and pp —?> $++$“ —>• i'^ip£~h's [111- Such an analysis was 

done in the context of the minimal type H seesaw model and the singly and doubly charged scalar bosons were taken 
to be mass degenerate. 

In the following analysis, we will consider that there is a charged scalar boson with a mass rrii above 300 GeV, 
which is consistent with measurement of the ATLAS and GMS collaborations. Since the form factors depend on the 
splitting between the masses of the charged scalar bosons Arriij = mj — rrii , we will use the parameters rrii and Am^ 
in our analysis below, along with the magnitude of the four-momentum of the virtual gauge boson, which we denote 
generically as ||p||. The region of interest corresponds to Arriij > 0 but for completeness we will also analyze the 
region —rrii < Arriij A 0 , namely, the scenario with rrij < rrii. Such a region, which corresponds to a very light 
charged scalar boson, is not favored by experimental data, but we will consider it in our analysis in order to show the 
consistency of our results. We will thus analyze the behavior of the form factors as functions of | \p\\ and Arriij for 
three illustrative values of rrii. We will show that all the form factors can have both real and imaginary parts. The 
latter appears when the magnitude of ||p|| reaches the value of the sum of the masses of the charged scalar bosons to 
which the virtual gauge boson is attached and it is a reflect of the fact that a pair of real charged scalars could be 
produced at a collider, rather than two virtual ones, via an off-shell gauge boson. 


A. The form factor /if 

To begin with, we show in Figure [5] the real (top plots) and imaginary (bottom plots) parts of the hf form factor 
as functions of the four-momentum magnitude ||p 21 | (left plots) and the splitting of the charged scalar boson masses 
Arriij (right plots). For best appreciation of the curves we show the absolute value of the real and imaginary parts of 
the form factor. We use three distinct values of rrii: 300, 400 and 500 GeV. 

We will first analyze the behavior of the real part of the /if form factor (top plots). It can be observed in Eq. 
([5]) that /if includes the term — p\ in the denominator and thus it becomes undefined when p\ —>■ which 
explains the sharp peak at ||p 2 || = ttiz in the top-left plot, where we include a vertical line for illustrative purposes. 
This effect is not in conflict with Landau-Yang’s theorem since the full vertex function ([3]) does vanish when all the 
Z gauge bosons are taken on-shell. In the ||p 2 || > 1600 GeV region of the top-left plot a dip appears in each curve, 
though in the plot it is only visible in the rrii = 300 GeV curve. This is due to a flip of sign of /if: in the ||p 2 || < rnz 
region the /if form factor is positive, whereas in the ||p 2 jj > 'mz region it is negative and changes sign again in the 
dip located at ||p 2 || > 1600. It is also interesting to note that the largest values of the real part of /if are reached 
around ||p 2 || mz. For instance, in units of lv[\{gfjgj^ ), we have |/if | ~ 10“^ for ||p 2 || — 100 GeV and |/if | ^ 10“® 
for IIP 2 II = 1000 GeV. This is to be compared with the average value obtained from the LEP lower and upper bounds 
on the /if form factor [l| (horizontal line of the top-left plot). As far as the behavior of /if as a function of Arriij is 
concerned, while the region with large Arriij corresponds to a heavy rrij, the —rrii A Arriij A 0 region corresponds to 
the scenario with a very light charged scalar boson with a mass lying in the interval 0 < rrij < rrii. Although this 
scenario seems to be ruled out by experimental data, we include it in our analysis for completeness. We note that 
in these and subsequent plots the stop of the curves is due to the reach of rrij = 0. As expected, since degenerate 
charged scalar bosons do not give rise to GP-violating form factors, the form factors vanish when Arriij = 0, this is 
why in the top-right plot we observe a sharp dip at Arriij = 0, which is due to the vanishing of the form factor. We 
also note that the real part of the /if form factor is not sensitive to a change in the splitting Arriij, so the different 
curves are almost indistinguishable. 

As far as the imaginary part of /if is concerned, which we show in the bottom plots of Fig. [SJ since the photon 








IIP 2 ll [GeV] 


Amij [GeVJ 


FIG. 5: Behavior of the real (top plots) and imaginary (bottom plots) parts of the hf form factor as functions of the four- 
momentum magnitude ||p 2 || of the virtual Z gauge boson (left plots) and the splitting of the charged scalar boson masses 
Arriij = rrij — rm (right plots). We use the indicated values of Arruj and ||p 2 || and each curve correspond to a distinct value 
of the charged scalar mass rrii-. rm — 300 GeV (solid line), 400 GeV (dashed line), and 500 GeV (dotted line). The horizontal 
line in the top-left plot corresponds to the average value of the LEP lower and upper limits [I[. We note that the curves stop 
because rrij = 0 is reached. 


must couple to the same charged scalar boson circulating into the loop, the Z gauge boson must necessarily couple 
to distinct charged scalar bosons. Thus the imaginary part of hf can only appear when ||p 2 || > Wi -|- rrij, which 
is evident in the curves of the bottom-left plot. It means that a higher energy would be required to measure such 
imaginary part unless there was a relatively light charged scalar, which is a scenario ruled out by experimental data. 
On the other hand, when the value of | |p 2 11 is fixed the imaginary part of hf can only appear in the small region 
0 < rrij < IIP 2 II ~ 'nii or —rrii < Arriij < IIP 2 II ~ as shown in the bottom-right plot. This interval becomes 
narrower for increasing rrii: for instance when ||p 2 || = 600 GeV, the imaginary part of hf is nonvanishing in the 
interval 0 < rrij < 300 GeV for rrii = 300 GeV, 0 < rrij < 200 GeV for rrii = 400 GeV and 0 < rrij < 100 GeV for 
rrii = 500 GeV. Again, we only include these results in our analysis for completeness. In general terms, we observe 
that the imaginary part of the hf form factor can have a size of similar order of magnitude than its real part, in 
the same interval of the region of parameters where the former is nonvanishing. However the maximal size of the 
imaginary part is reached for a very light charged scalar, whereas the maximal size of the real part can be reached 
around the Z resonance, where there is no imaginary part. 


B. The form factor hf 

We now present in Figure [S] the corresponding plots for the hf form factor, namely, we show the behavior of the 
real (top plots) and imaginary (bottom plots) parts of the form factor hf as functions of | |p 211 (left plots) and Arriij 
(right plots). We use the same set of parameters as in Fig. [SI We note that in general the hf form factor shows 
a similar behavior to that of the hf form factor, though there are some slight differences. We first analyze the top 
plots of Fig. ini which show the real part of hf. We observe that apart from the sharp peak at ||p 2 || = "mz in the 
top-left plot, there are also dips at high energy, such as occurs in the respective hf curves. Such dips, which are a 
result of the flip of sign of the form factor, are now shifted to the left and they appear at a higher ||p 2 || for a smaller 
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rrii, which is opposite to the behavior of hf. Thus the dip for rrii = 300 GeV curve is not shown in the plot. In 
the IIP 2 II < fnz region /if is negative, whereas in the IIP 2 II > "m-z region it is positive, contrary to what happens 
with /if. After the dip at ||p 2 || > 1000 GeV, /if becomes negative again. We also note that the magnitude of /if is 
greater than that of /if, although for very high values of ||p 2 || or very heavy rrii, the both /if and /if are considerably 
suppressed. As expected, the largest values of the real part of /if are reached around ||p 2 || mz- For instance, in 
units of we have |/if | ~ 10“® for ||p 2 || = 100 GeV and |/if | ~ 10“^ for ||p 2 || = 1000 GeV. It can also be 

observed that, /if is more sensitive than /if to a change in the value of rrii. Regarding the top-right plot, again we 
observe that the /if form factor vanishes when rrii = frij, which shows the consistency of our result. The curves of the 
top-right plot show a dip in the —rrii < Arriij < 0 region , which are due to a change of sign and do not appear in the 
case of the /if form factor. Although /if can reach its largest values in this region, as explained above, it corresponds 
to the case of a very light charged scalar. 


WpiW [GeV] 


Aniij [GeVJ 



FIG. 6: The same as in Fig. [5] but for the /if form factor. 

As in the case of /if (both hf and /if form factors arise from the same vertex function), the imaginary part 
of /if would be nonvanishing when ||p2|| > rm + rrij, which again is evident in the bottom-left plot of Figure |6l 
Therefore, to observe such an imaginary part a higher energy than that required to observe the respective real part 
would be required. Furthermore, although the sizes of the imaginary part and the real part are similar, the real part 
can reach larger values in the interval where the imaginary part vanishes. Thus, in general the imaginary part is 
smaller than the largest possible values of the real part, which is reached around the Z resonance. As far as the 
behavior of the imaginary part of /if as a function of Arriij is concerned, as expected it is nonvanishing in the interval 
—rrii < Arriij 'A IIP2II ~ 2mi, where rrij is very light. In this region the imaginary part of /if can reach values slightly 
larger than the real part. 


C. The form factor /f 

We now turn to analyze the fj form factor. We note that it becomes undefined when Hgjl —>■ 0 as can be inferred 
from Eq. m- Nevertheless, Landau-Yang’s theorem is obeyed as can be deduced after a closer inspection of Eq. 
dH). In the top plots of Figure[7]we show the real part of f] as a function of the four-momentum magnitude ||g|| (left 
plot) and the mass splitting Arriij (right plot). On the other hand, similar plots are shown in the bottom part of the 
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panel that illustrate the behavior of the imaginary part of f2- In this analysis we will use the same set of values 
used in Fig. [5] 


llqll [GeV] 


Aniij [GeV] 



llqll [GeV] 


Aniij [GeV] 


FIG. 7: Behavior of the real (top plots) and imaginary (bottom plots) parts of the fj form factor as a function of the four- 
momentum magnitude ||?||of the virtual photon (left plots) and the splitting of the charged scalar boson masses Arriij (right 
plots). We use the indicated values for Arruj and ||g|| and each curve correspond to a distinct value of the charged scalar mass 
rrii'. rrii = 300 GeV (solid line), 400 GeV (dashed line), and 500 GeV (dotted line). As explained in the text, the imaginary part 
of f2 arises only when ||ij|| > 2mi in the bottom-left plot and 0 < rrij < or —rrii < Arriij < — rrii in the bottom-right 

plot. 


As far as the real part of the f2 form factor is concerned, we note that the dip appearing in each curve of the 
top-left plot corresponds to a flip of sign of f2, which turns from negative to positive at ||( 7 || ~ 950, 1100, and 1500 
GeV for mi =300, 400, and 500 GeV, respectively. We also observe that this form factor decouples at high energy, 
where it has a negligible magnitude, but in the interval between 100 GeV and 900 GeV it can reach values as higher 
as 10“® — 10“®, in units of ). It is also interesting to note that the real part of f2 is not very sensitive to a 

variation of the charged scalar mass m^, which contrasts with the behavior of the imaginary part. On the other hand, 
in the top-right plot of FigureO we show the behavior of f2 as a function of for ||g|| = 600 GeV and the same 
three values of used in previous analyses. We observe that the magnitude of f2 does not increase significantly as 
Amij increases. As expected, this form factor vanishes when Amij = 0. We note that the maximal values of the real 
part of f2 are of the order of 10 “® — 10 “^ times the factor lin{gfjg2*). 

In the bottom plots of Figure [7] we show the behavior of the imaginary part of f2- Since the virtual photon 
must necessarily couple to the same charged scalar boson, the imaginary part of this form factor would arise when 
IIqII > min (2mi, 2mjj), which is evident in the bottom-right plot, where the imaginary part is nonzero for ||q|| > 2 mi. 
For fixed m^, the interval for nonvanishing imaginary part can be written as —mi < Amij < —or 0 < mj < 

This is the reason why the curves in the bottom-right plot are nonvanishing only in the small region where —300 
GeV < Amy < 0 GeV for = 300 GeV, -400 GeV < Amy < -100 GeV for m, = 400 GeV and -500 GeV 
< Amy < —200 GeV for m^ = 500 GeV. Again we note that the imaginary part of f2 has a similar order of 
magnitude than its real part, namely 10“® times Imlgfjg^i*). However, as in the previous form factors, the real part 
can reach larger values in the region where the imaginary part vanishes than in the region where both of them are 


nonzero. 
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D. The form factor // 

Finally we show in Figure [8] the real (top plots) and imaginary (bottom plots) parts of the form factor as a 
function of the four-momentum magnitude HgH (left plots) and the mass splitting Arriij (right plots). We consider 
the same scenarios analyzed above. 


I [GeV] 


Amij [GeVJ 





I" 


1 [GeV] 


Amij [GeV] 


FIG. 8: The same as in Fig. [7|but for the ff form factor. Unlike the f2 form factor, which becomes undefined : 
/f form factor gets undefined in ||g|| = mz- 


= 0, the 


We first analyze the behavior of the real part of ff (top plots) as a function of ||g|| and Arriij. This form factor has 
a sharp peak at ||(?|| = mz due to the factor — m\ in the denominator. There are also two dips appearing in the 
curves of the top-left plot, which as in the previous cases are due to a sign flip of this form factor. We observe that, 
unlike the /if 2 form factors, // flips sign twice for mz < ||g|| < 2000 GeV. One of such sign flips occurs at H^H ~ 110 
GeV, regardless of the value, and the second flip occurs at HgH ~ 800, 1000, and 1250 GeV for = 300, 400, and 
500 GeV, respectively. Inside the region enclosed by the two dips the real part of is positive, and it is negative 
outside this region. Although the larger values of fi are reached around the Z resonance, there is an increase of the 
real part of in the interval between ||g|| =140 GeV and ||g|| =700 GeV, where the real part of goes from 10“® 
up to I 0 -^ in units of Im(g,f g?*). 

The /f for factor develops an imaginary part when | |q| | > min (2mi, mi + mj, 2mj) since in this loop the virtual Z 
boson can now couple to the same or distinct charged scalar bosons. In the bottom-left plot we have m_,- > mi, thus 
the imaginary part of ff arises when ||q|| > 2mi. On the other hand, for the chosen values of mi, in the left plot we 
have mi > thus the interval where the imaginary part of ff is nonzero is given by ||( 7 || > 2mj or in terms of the 

mass-splitting —mi < Amij A — Wi, as observed in the curves shown in the bottom-right plot. Thus fj and ff 
develop imaginary parts in the same interval. Most part of this interval correspond to a very light charged scalar and 
as occurs with the other form factors, the imaginary part of ff has a size of the same order of magnitude than its 
real part, which however can reach much higher values outside this region. 

It is interesting to note that the // form factor was also studied in the framework of a THDM where the respective 
contribution is induced by three nondegenerate neutral scalar bosons [ 2 ^. In such a model three different nondiagonal 
complex couplings On arise in the neutral scalar sector. It was reported in 221 that ff can reach values from 
10“® X O11O12O13 to I0“® X O11O12O13, where the following set of values for the free parameters was used: = 
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llgll = 200 GeV, M 3 = 250 GeV, M 2 = 150 GeV, and 60 GeV < Mi < 150 GeV, with Mi the masses of the neutral 
scalar bosons. In this case the O 11 O 12 O 13 factor could suppress considerably the THDM contribution to //, just as 
happens with the contribution of our charged scalar bosons, which could be suppressed by the factor ). 


V. CONCLUSIONS 

We have presented an analysis of the one-loop contributions to the GP-violating form factors associated with the 
ZZ*^, ZZ^* and ZZZ* couplings in the framework of an arbitrary effective model with at least two nondegenerate 
charged scalar bosons that couple nondiagonally to the Z gauge boson. Such form factors are induced only when the 
nondiagonal coupling constants gfj are complex and have an imaginary phase. Our analysis is independent of any 
specific value of the coupling constants, so our results are scaled by the coefficient (/him((/^^(/?*). We considered a 
charged scalar boson with mass above 300 GeV, which is consistent with experimental constraints, and analyze the 
behavior of the real and imaginary parts of the form factors as functions of the four-momentum of the virtual gauge 
boson I Ip 11 or ||g|| and the splitting of the masses of the charged scalar bosons Annj = mj — iru. Although the region 
which is consistent with experimental data corresponds to Arriij > 0 , for completeness we also include the results for 
—TTii < Arriij < 0, which corresponds to nij lighter than rm . As far as the orders of magnitude of the form factors are 
concerned, they are as follows in units of (/him((/^^(/?*): |hf | ^ 10 “^ — 10 “^, |hf | ~ 10 “^ — 10 “®, I /4 | ~ 10 “® — 10 “^ 
and I//I ~ 10“® — I0“®, for values of ||p|| about a few hundreds of GeVs. When the magnitude of ||p|| increases, the 
real part of the form factors decreases smoothly and decouple at high energy. Furthermore, it is not very sensitive to 
large values of Arriij. As for its imaginary part, for an intermediate value of ||p|| it can arise in a small region of values 
of rrii and rrij, which corresponds mainly to the scenario with a relatively light charged scalar. Such scenario is not 
favored by current constraints on the mass of a charged scalar boson. On the other hand, for heavy nii and rrij, the 
imaginary part of the form factors can only arise for large ||p||. Therefore a higher energy would be required in order 
that the form factors could develop an imaginary part, which can have a magnitude similar to the corresponding real 
part, though the latter can reach larger values in the region where the former vanishes. We also find that, except for 
a proportionality factor, our results for the ff form factor are of the same order of magnitude than the contributions 
arising in a THDM with three nondegenerate neutral scalar bosons that couple nondiagonally to the Z gauge boson. 
It is worth noticing that our results are valid for the contribution of nondegenerate doubly charged scalar bosons, 
though in this case the form factors associated with the ZZ*j and ZZj* couplings would get an additional factor of 
two. 

Although the LHC was down temporarily, it has resumed their operations with a higher center of mass energy 
and a higher integrated luminosity. A significant improvement in the experimental limits of the TNGBCs is thus 
expected in the forthcoming years. Therefore it is necessary to examine any potential contribution to the respective 
form factors. 
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